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ACOUSTIC SCATTERING AND ABSORPTION BY A RIGID 
POROUS ELLIPTIC CYLINDRICAL SHELL* 

By William E. Zormnski 
Langley Research Center 

SUMMARY 

Solutions for scattering of a plane wave from a thin rigid porous elliptic cylindrical 
shell are found in terms of Mathieu fimctions. The scattering problem reduces to an infi- 
nite set of linear equations for the coefficients of the velocity expansion on the cylindrical 
surface. Coupling terms in these equations are given by integrals which are evaluated by 
exact methods. Solutions are found by the method of reduction for the internal standing 
waves and the external scattered waves, A special case of variable acoustic impedance 
is foxmd where the infinite systems uncouple so that closed-form solutions may be obtained 
for the velocity coefficients. Computations are made for scattering in the intermediate 
frequency range, where the wavelength is of the same order of magnitude as the major 
axis of the cylinder. Tabulations of the scattered and dissipated energies are given for 
various cylinder eccentricities, impedances, angles of incidence, and frequencies. 

INTRODUCTION 

The problem of acoustic scattering from a porous cylinder originates from the con- 
temporary need for methods to reduce noise from aircraft turbofan engines. In refer- 
ence 1 Marsh has shown that a practical method for reducing noise from turbofan engines 
is to install "broad-band resonators" inside the engine nacelle. Because of considerations 
of weight, safety, and endurance, these resonators are usually made of thin porous sheets 
of material (either metallic or fiber -glass — plastic) which are fastened to a honeycomb 
wall structure. The cavities behind the porous sheet are usually about one-quarter wave- 
length deep, since this depth gives good absorbing qualities. 

In general, the greater the exposed area of porous material, the more the soimd is 
absorbed, so that engine designers must look for ways to alter the engine geometry to 
increase this area. Of course, this increase must be accomplished without upsetting the 
basic flow field within the engine, which presumably has already been optimized on a 

*The material presented herein is based on a thesis entitled "Acoustic Scattering by 
a Porous Elliptic Cylinder With Nonlinear Resistance" submitted in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy in Engineering Mechanics, 
Virginia Pol 5 d;echnic Institute, Blacksbtirg, Virginia, March 1970. 



performance basis. One approach to this problem has been to install thin porous double - 
walled radial spokes or circumferential rings in the engine. The optimum size (from the 
acoustical viewpoint) of these devices is not known, but it is reasonable to assume that 
their dimensions are of the same order of magnitude as the wavelength of the dominant 
tone. In cross section, such a spoke or ring would appear roughly as an ellipse. 

The literature on diffraction is extensive. The reader interested in this subject 
should consult the important review article by Bouwkamp (ref. 2). More than 500 papers 
published between 1940 and 1954 are reviewed in Bouwkamp 's article, and the application 
of Mathieu functions to diffraction problems is discussed. Briefly, Mathieu functions may 
be used in the study of diffraction (or scattering) from elliptic cylinders. The article 
indicates that research prior to 1954 was limited to studies of the problem of scattering 
from a strip or slit (cylinder with eccentricity 1). This scattering problem has also been 
studied with an integral equation formulation (ref. 3). 

Since Bouwkamp's survey, several papers have been published which deal with cyl- 
inders having eccentricity less than 1. In 1963 Yeh discussed the problem of a penetrable 
strip in terms of a Mathieu function series (ref. 4), and Barakat (ref. 5) made a study of 
the elliptic cylinder with various eccentricities. Barakat considered the exterior problem 
by using the classical Dirichlet (sotmd soft) and Neumann (sound hard) botmdary conditions. 
This work made use of recently computed tables by Barakat, Houston, and Levin (ref. 6). 
These tables are appropriate for scattering problems in the long -wavelength range. In 
1964 Burke and Twersky (see refs. 7, 8, and 9) began investigations in which low-frequency 
approximations were used. 

In reference 7 Burke considered the problem of a cylinder composed of fluid with a 
density different from that of the medium in which it is immersed. Boundary conditions 
at the surface of the elliptic cylinder were continuity of normal velocity and continuity of 
pressure. His paper, like Yeh's (ref. 4), considered the coupled interior and exterior 
cylinder problem. The analysis of this problem resulted in an infinite set of algebraic 
equations for the scattered waves. Burke indicated the connection between truncating 
these sets of equations and long -wavelength approximations. His series solution gave the 
far -field scattering amplitude up to the sixth power of the wave number, and the near -field 
amplitude (internal and external) up to the third power. Terms in the truncated equations 
could be evaluated from the previously mentioned tables of Barakat, Houston, and Levin 
(ref. 6), from the tables of Wiltse and King (refs. 10 and 11), or from computations based 
on known expansions of the Mathieu functions in terms of the other tabulated transcen- 
dental fimctions. 

This paper considers the scattering of a plane wave from a porous cylindrical shell. 
The cylinder is a thin elliptic shell which is made of a rigid porous material. The solu- 
tions both inside and outside of the cylinder are obtained in a general form and expressed 
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in terms of the unknown velocity normal to the cylinder surface. It is assumed that this 
normal velocity is continuous at the cylinder surface and that the pressure differential is 
related to the normal velocity by the complex impedance of the porous material. This 
relation between pressure differential and velocity is used to find the solution for the 
velocity. With the velocity through the cylinder surface known, calculations may be made 
for other quantities such as scattered energy and pressures inside and outside of the 
cylinder. 


SYMBOLS 


Ar,Br,Cr,Dr coefficients 

b,bej.,bOr characteristic values 

C^(x) Gegenbauer polynomial 

c speed of soimd 

De^^^ Fotirier coefficient of even periodic Mathieu fimction 

Do^ ’ Fourier coefficient of odd periodic Mathieu fimction 

DR ratio of dissipated to incident power 

e base of natural system of logarithms, 2.718 

g(^,rj) = h/cosh^l - cos^ry 
Her^^\s,|) = Jer(s,|) + iNer(s,?) 

Hej.^^\s,|) = Jer(s,?) - iNer(s,4) 

HOr^^\s,|) = JOj.(s,^) +iNoj.(s,^) 

Hor^^^s,^) = JOr(s,4) - iNor(s,|) 
h dimensionless focal length 

i unit imaginary number 
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Jer(s,|), JOr(s,^) radial Mathieu functions of the first kind 

JnjYn Bessel functions of integer order 

j,k,Z,m,n,p,r,v integers 

Se^(s,??)Ser(s,r7)dr? 

= |e(^o»^)l ^So^(s,??)SOr(s,r7)dr7 

p2ir p -,2 

Nr = j^ [_Ser(s,?7)J d?7 

r>27T p -i2 

Nr = J^ lSor(s,r7)J d?7 

Nej.(s,^),Nop(s, radial Mathieu functions of the second kind 
P dimensionless power 

p dimensionless pressure 

q dimensionless velocity 

component of velocity normal to surface, | = Constant 
real part of a complex number 


Re( ) 


Ser(s,rj) even periodic Mathieu function 


Sor(s,? 7 ) odd periodic Mathieu function 


SR 


ratio of scattered to incident power 


s = h" 


dimensionless time 


Uj.,Vj. coefficients in series for 
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w projected width of cylinder 

x,y dimensionless rectangular coordinates 

Z specific impedance of cylinder 

a angle of incidence 

7} circumferential elliptical coordinate 

^ radial elliptical coordinate 

coordinate of cylinder surface 
p mass density 

(p,xj^ solutions to Mathieu's equation and Mathieu's modified equation, respectively 

(t) circular frequency 

Subscripts: 

d dissipation in porous material 

1 incident 

p incident plane waves 

s scattered waves 

Superscripts: 

i inside of cylinder 

o outside of cylinder 

(1) denotes diverging waves 

(2) denotes converging waves 
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complex conjugate 


* 


A single prime denotes a first derivative (except for Np, which has been defined 
separately). Double primes denote a second derivative. 

ANALYSIS 


Governing Equations 

Figure 1 shows the elliptical cylinder coordinates (| and 77 ) which will be used to 
study scattering from an elliptic cylinder. These coordinates are defined by 


X = h cosh I cos 7} 


y = h sinh ^ sin rj 


( 1 ) 


As shown in figure 1, curves with ^ equal to a constant are ellipses with focal points 
at X = ±h. Curves with rj equal to a constant are hyperbolas with focal points which 
correspond to those of the ellipses. 



Figure 1.- Elliptical coordinates. 
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The problem being considered herein is the scattering of a plane wave from a 
porous cylinder whose surface corresponds to the elliptical coordinate This cyl- 

inder has a semimajor axis h cosh and semiminor axis h sinh The elliptical 
coordinates are naturally dimensionless quantities, and it is convenient to introduce 
dimensionless parameters for the other variables. Thus, h, p, q, t, x, y, and Z 
are all defined as dimensionless parameters: 


h = 


O) X Focal length 
c 


P = 

q = 


Pressure 

pc2 

Velocity 

c 


and 


t = CO X Time 

CO X Longitudinal coordinate 


X = 


y = 


z = 


CO X Lateral coordinate 


Impedance 

pc 


The acoustic field equations are 

V^p + p = 0 


q = -iVp 


( 2 ) 

(3) 


where the time factor exp(-it) is understood. 

The pressure and its gradient must be continuous inside and outside of the cylinder. 
Therefore, there are two continuity conditions at | = 0: 


pM 0,?7) = p^(0,-?7) (4a) 

^pH0,rj) ^ -8pi(0,-T?) , . 

At the surface of the porous cylinder, where | there are two boundary conditions 

which relate the internal and external pressure fields. In order for mass to be conserved, 
the acoustic velocity normal to the cylinder sm-face must be continuous; that is. 
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( 5 ) 


However, there will be a pressure drop at the cylinder surface due to viscous dissipation 
and inertial effects within the porous material. This pressure drop is related to the 
acoustic velocity through the specific impedance of the material Z; that is, 

The field equations in elliptical coordinates become 

+ h^(cosh^| - cos^?7)p = 0 (7) 


and 

= -i[g(«,T7)] H 

Separating variables in equation (7) by letting 


gives 


4>"{q) + ^b - s cos^'r^4>h]) ~ 0 


( 8 ) 


( 9 ) 


(10) 


and 


- (b - s cosh^^)i^(^) = 0 


( 11 ) 


where 


s = h^ (12) 

Since the dimensionless parameter h is proportional to frequency, s is proportional 
to frequency-squared. 

Equations (10) and (11) are called Mathieu's equation and Mathieu's modified equa- 
tion, respectively. The solution <p of equation (10) must be a periodic function of q, 
of period tt or 2n; the parameter b, which is a function of s, must be one of a coxint- 
ably infinite set of characteristic values for every s. Fourier series representations 
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of these periodic solutions are given by equations (40) to (43) in the appendix. These 
solutions form a complete orthogonal set 



Seris,v)Sei(s,7j)dr] 


0 





(^ = r) 


(13a) 



Sor(s,77)So^(s,r7)d?7 = 0 



? (13b) 

(l = r)J 


For each characteristic value, there are two independent solutions of equation (11). These 
solutions are called radial solutions, or hyperbolic cylinder functions. Series represen- 
tations of these solutions are given in reference 12. 


Solution for Scattering From Cylinder 

In terms of the notation of reference 12, the general equation for the pressure field 
is 


r=0 

[b 




Ser(s, 7j) 




Soj.(s,?7) 


(14) 


The asymptotic expansions of the radial solutions show that solutions with the 
superscript (1) are diverging waves while the superscript (2) indicates converging waves. 

It is often convenient to work with the real and imaginary parts of the radial solu- 
tions. In reference 12, the radial solutions are given as 

Her^^^(s,|) = Jer(s,^) + iNer(s,?) (15a) 


and 



(s,5) = Jor(s,|) + iNor(s,|) 


(15b) 
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I iiiiiniii 


III 


The Wronskian of these solutions is 1; that is 


JOj'Nej’ “ Nej’Jeji' — 1 


and 


JOj.Noj. - NOj.JOj. = 1 


(16a) 


(16b) 


It can be shown that, in order for the pressure and velocity to be continuous at | = 0, the 
pressure inside the cylinder must be 


pH^,v) = ^ jArJer(s,^)Ser(s,?7) + BrJor(s,|)Sor(s,r7) 
r=0 


(17) 


From succeeding equations, the dependence of the functions on the parameter s can be 
understood. 

It is convenient to formulate the problem in terms of the normal velocity throi^h 
the cylinder surface. An expansion for this velocity may be taken in the form 

OO 

^ |UrSer(??) + Vj.Soj.(r?^ (18) 

r=0 


so that a solution for the coefficients Uj. and Vj. is equivalent to finding the velocity 
distribution. With equations (8) and (18), it may be shown that equation (17) becomes 


u>. \ -V Hr'l®r(?o) c, f ^ Vj.Jor(^o) c, , ^ 

^ / ~r"Tt \ \ ■■■ Sor(??) 

r=oL'^®r(^o) J°r(y 


(19) 


The pressure due to plane waves of unit amplitude moving along a line which is 
inclined at an angle a to the major axis may be represented by (refs. 13 and 14) 


Pp(^>v) = exp[i(x cos 0 ! + y sin a;)] = ^ ^ (i)^ 


r=0 


Se„(Q!) . , So,,(q!) 

JerU)Sej,(r]) + — ^ — Joj.(^)SOj.(?7) 
iN V 


Nr 


(20) 


where (x cos O! + y sin a) is the normal distance from the origin to the wave front. 
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The pressure due to waves scattered from the cylinder may be given as 


Pg(|,77)= ^ [crHer^^\?)Ser(r?) + DrHor^^\^)Sor(T?)] 


( 21 ) 


r=0 


The total pressure outside of the cylinder is the sum of equations (20) and (21); that is, 

Pr> + Po* If equation (8), used to calculate the normal velocity from the total pressure, is 
matched with the velocity expansion (eq. (18)), Cj. and Dj. may be eliminated from the 
expression for the external pressure, which becomes 




r=0 


He, 


( 1 ) 


He;(l) 


(^o) 


l^(i)^ 


(^q) 

Sej.(a!) 


Sep (77) + Vr 


Ho, 


( 1 ) 


(^o) 


«Or‘**({o) 


Sor(rj) 


NrHer(l)(lo) 


, , SOt.(0!) , , 

Ser(^) + Sor(r?) 


n;ho;<«(4o) 


( 22 ) 


The interior -pressure solution (eq. (19)) and the exterior -pressure solution (eq. (22)) 
are each expressed in terms of the velocity at the cylinder surface. The pressure differ- 
ential and velocity at the surface must satisfy equation (6) so that substituting equa- 
tions (19), (22), and (18) into equation (6) and utilizing the orthogonality properties of the 
periodic Mathieu functions (Sej. and Soj.) give the four infinite sets of equations 


Z ^ Moj ‘j-K.Uov. + 


r=0 


‘^2l,2r^2r 


% _ i^^'*'^l/^Se2^(o!) 

ill ~ /1\ 

■Iea(io)Hek' '(lo) 


(23) 


and 


N2i+lU2i+l j2^+2|/^Se2^H-l(«) 


r=0 


r=0 


^ ^ %+2,2r+2^2r+2 + — 


r=0 




Il2i+lV2i+l 

i2Z+2|/3^So2^^j(o!) 

J°k+l(^o)II°k+l^^^(^o) 


'^21+2^21+2 

i2^+^|/§^So2^+2(") 

J°a+2(So)Hoa+2‘^’(5o) 

Hok+2‘‘’(«o) 


(24) 


(25) 


(26) 
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where I = 0,1,2, . . . and where and are integrals that are defined in the 

symbol list and evaluated in the appendix. These four infinite sets of equations give the 
coefficients of the velocity expansion on the cylinder surface (Up and Vp) . Once the 
velocity at the cylinder surface is obtained, the pressure or velocity at any point may be 
computed. 


A Special Case of Variable Impedance 

There is a special case of variable material impedance for which it is possible to 
obtain a simple algebraic expression for Up and Vp instead of the infinite sets of 
equations. If 


Z = 


where Z' is a constant, then 


and 


Up = 


F+lfeep(o!)Je;(^o) 
Np[l + Z’Je;(|o)He;.(l)(lo) 


Vp = 


tr+l 


l^SOp(Q!)Jor(4o) 


n; 


1 + Z’Jo’ 


(27) 


(28) 


(29) 


Scattering Formulas 

Some quantities of interest are the incident power, the scattered and dissipated 
power, and the far -field root-mean-square pressure. The derivation of formulas for 
these quantities is tedious, but routine, and only the definitions and final results are given 
herein. 

The projected width of the cylinder on the oncoming plane wave is 

w = (30) 


as shown in figure 1. Consider a section of the cylinder whose dimensionless length 
X Lengtl^ power due to the incident wave which passes through the projected 

area of this unit length of the cylinder is 
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p _ w j?c dt - ^ 


(31) 


In equation (31) a xinit dimensionless length is being considered, 
this length is found by multiplying equation (31) by the factor ^ 

IS 

Pg = ^ Re 

It may be shown that equation (32) reduces to 

OO 


r=0 


Cr ^Nr + iDrl^Nr 


where 




Hei.(l)(^o) 


iUr - iri/^5^Je;(?o) 


and 


Dr = 






The power dissipated in the porous material is 

■^2 tt 




z 


jUrlX 


IVrl^N^. 


r=0 l^®r(^o)^®r^ \^o) '^°r(^o)^°r^ ^(^o 


(1). 


+ i^+l|/87 


Ser(Q!)Ur ^ Sor(o!)V^ 




The actual power for 
5 

The scattered power 
(32) 


(33) 


(34) 


(35) 


} 

(36) 


13 



At large distances from the cylinder where ^ » 1, the rms pressure is given by 


ihe^] 


P 


rms 


WJ 

^ (~i)^ |pr®®r(^) “ ®r+l®*^r+l(^^ 


r=0 


(37) 


RESULTS AND DISCUSSION 

Scattering computations have been made for wavelengths that are of the same order 
of magnitude as the major axis of the cylinder. The longest wavelength considered is 
about three times the length of the major axis. The parameter s, which is proportional 
to the square of the frequency, is 4 for this low-frequency case. The shortest wavelength 
is about one -third the length of the major axis, which corresponds to s = 100. 

The parameter determines the eccentricity of the cylinder. The ratio of the 
minor axis to the major axis of the cylinder is approximately equal to when 
small. Computations have been made for values of 0.4, and 0.8. The largest 

value of gives a cylinder with a minor- to major -axis ratio of about 0.68. 

For the frequencies considered in this paper, the Mathieu fimction series for the 
pressure and velocity fields were always sufficiently accurate with 14 terms, so it was 
necessary to evaluate Mathieu functions from order 0 to order 14. 

Evaluation of both the periodic and radial Mathieu fimctions (ref. 12) depends on 
computing the coefficients De^^^ and E)o^^\ which appear in the Fourier series given 
in the appendix. These coefficients are functions of the parameter s. 

Evaluation of Mathieu Functions 

Although the coefficients for the Mathieu functions are tabulated in reference 12 for 
the range of s used in this paper, it was easier to program the computer to use the 
characteristic values bej. and boj. to compute these coefficients rather than to read 
and interpolate the tables of coefficients. The method given by Blanch in reference 15 
was used for this computation. Results were checked against the tables of coefficients 
and were found to agree to eight significant figures. With this method of computation it 
is also possible to obtain the coefficients for s > 100 (which are not tabulated) from the 
tables of characteristic values for large s (s > 100) in reference 12. With the coeffi- 
cients available, Fourier series could be used for computation of the periodic solutions. 

The evaluation of the radial Mathieu fimctions depends on computing and Yn 
(Bessel functions of the first kind and second kind, respectively). The Bessel functions 
of the first kind were computed by Abramowitz’s method (see ref. 16). In this method a 
unit value is assigned to Jjj (where k is large compared with the argument of the 
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Bessel fmiction), is taken as zero, and then backward recurrence is used to find 

Jjj for n < k. The normalization condition is then used to adjust the magnitude of the 
Jn's. The Bessel fimction Yq is obtained from its series representation in terms of 
Jji, Yj is fo\md from the Wronskian relation, and the rest of the Y^'s are found by 
forward recimrence. Bessel fimctions computed by this method were accurate to nine or 
ten significant figures for the range of argmnents used. 

Because of the scarcity of tabulated data available on the radial Mathieu functions, 
these fimctions were calculated by two different methods. The Bessel fimction product 
series (eqs. 3.05 to 3.13 in ref. 12) was used; in addition, the Bessel function series 
(eqs. 3.03, 3.04, 3.15, and 3.16 in ref. 12) was used for comparison. The latter series 
involves Bessel functions with argument i/s' cosh Each series was checked at the 
origin = oj against the joining factors which are tabulated in reference 12. Also, each 
series was checked against the tables of Barakat, Houston, and Levin (ref. 6) and against 
the tables of Blanch and Clemm (ref. 17) for s = 4. This value was the upper limit 
(on s) of the available tables and the lower limit on s for this study. 

The radial Mathieu functions of the first kind were the most difficult to evaluate. 

The Bessel function product series for Jeo(4,0) checked against the National Bureau of 
Standards tables (ref. 12) to seven significant figures. However, this accuracy decreased 
to one significant figure for Je0(4,O), and functions of higher order could not be computed 
by this method. The accuracy of the product series for the Jej.(s,^) and JOj.(s, 4) 
functions increased with increasing s; JeQ(100,0) was accurate to nine figures, and 
Jej 4 ( 100 , 0 ) retained a two -figure accuracy. 

The Bessel fimction series proved to be a superior method for computing JCj. and 
Joj.. By this method Jeo(4,0) was accurate to seven significant figures. Although this 
result was not an improvement over that obtained by using the product series, Je0(4,O) 
was found to be accurate to six significant figures, which represented an improvement of 
five significant figures, and Jeg(4,0) was accurate to five significant figures. However, 
Jeio(4,0) was accurate to only one significant figure, and Jei2(4,0) and Jei4(4,0) 
could not be computed. The accuracy of the Bessel function series for Jej.(s,0) also 
tended to improve with increasing s; Jeo(100,0) was accurate to six figures and 
Jei4(100,0) was accurate to eight figures. Note that the higher order functions are more 
accurate than the lower order functions when s is large. This turns out to be a very 
important property, since, in the scattering problem, the long -wavelength (s = 4) solutions 
require only the lower order terms for convergence. The short -wavelength solutions 
(s = 100) require more terms, but it is no problem to compute these terms for the large 
values of s. 

The radial Mathieu functions of the second kind were computed by the Bessel func- 
tion product series. For 0 = r £ 15, Nej;.(4,0) and Nor(4,0) were accurate to eight 





significant figures. The accuracy of the product series for Nej.(s,0) and Noj.(s,0) 
decreased with increasing s. For example, Neo( 100,0) was only accurate to three 
significant figures; Nei 4 ( 100 , 0 ) had six -figure accuracy, a loss of only two places. 

The accuracy of both methods of computation improved with distance from the ori- 
gin (increasing This fact was established by comparing the results of each series 

and assuming that, if they were the same, they were correct. Larger values of | 
correspond to smaller values of cylinder eccentricity; thus it is generally easier to eval- 
uate the radial fimctions on a nearly circular cylinder than on a nearly flat cylinder. 

In the problem of scattering from a porous cylinder, the Bessel function product 
series was used to evaluate the functions of the second kind and their derivatives. The 
Bessel function series was used to find the functions of the first kind, and the derivatives 
of these functions were found from the Wronskian relations (eqs. (16)). 

The magnitude and phase of the complex radial Mathieu functions Hej.(s,^) and 
Hoj.(s, ^) could be computed to eight significant figures for s = 4, and to six figures for 
s = 100. This fact is somewhat siirprising in view of the preceding discussion of the 
accuracy of the radial Mathieu functions of the first and second kinds which are, respec- 
tively, the real and imaginary parts of the complex function. The inaccurate (real) part 
of the complex number is much less in magnitude than the accurate (imaginary) part; thus, 
the complex representation of the number is fairly good throughout the range of s. 

The evaluation of the coefficients in the scattering equations (eqs. (23) to (26)) 
depends directly on the evaluation of the Mathieu functions. The coefficients Uj. and Vj. 
in equations (23) to (26) are only as accurate as the radial functions of the first kind (jep 
and Joj.) and their derivatives (Je^ and Jo^.). These quantities could be computed to at 
least five significant figures for s = 4 and orders up to eight. At s = 100 the fimctions 
of the first kind could be determined to at least six figures for all orders, but the func- 
tions of the second kind (Ncj. and Noj.) could be obtained to only three figures for the 
lower orders. Since the Wronskian relations were used to solve for the derivatives, the 
derivatives of Jej. and Joj. are limited to three -figure accuracy. The acciiracy prob- 
lem may be associated with the waves within the cylinder since this field is given in terms 
of radial functions of the first kind. 

The nonhomogeneous terms in equations (23) to (26) may be computed to six signifi- 
cant figures since they involve only the complex radial functions and the periodic func- 
tions. The normalizing factors N^. and are computed from the coefficients Del^^ 
and > and so are accurate to eight figures. 

Computation of the coupling coefficients M^ ^ and is discussed in the appen- 

dix. The series presented therein were used to compute these coefficients. The inte- 
grals l2m were computed to eight significant figures to obtain eight -figure accuracy for 
the coupling terms. 
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The four infinite sets of equations (eqs. (23) to (26)) were solved approximately by 
truncating the equations to a finite set. This method of solution, which is also called the 
method of reduction, is discussed in reference 18. Computations were made by using 
matrices of orders 5, 6, and 7, which gave 10, 12, and 14 terms, respectively, in the 
series for q^(4Q,T?)« K was observed that the numerical results did not change signifi- 
cantly as the order of the matrices increased to 7; this fact implies convergence of the 
solution by the method of reduction. 

Scattering From a Porous Cylinder 

Focal lengths h of 2, 4, 6, 8, and 10 were used to determine the effect of frequency 
on scattering and dissipation. These values correspond to lengths of the major axis of 
about one -half wavelength to over three wavelengths. Solutions were computed for real 
values of Z of 1, 2, and 4. The angle of incidence a was varied in 30° increments 
from 0° to 90°, and the radial coordinate was taken as 0.2, 0.4, and 0.8. Table I 
gives the ratios of scattered and dissipated power to incident power (SR and DR) for 
these conditions. 

TABLE I.- SCATTERED AND DISSIPATED POWER RATIOS 


(a) h = 2.0 



ot, 

SR 

DR 

SR 

DR 

SR 

DR 


deg 

for 2 

= 1.0 

for 2 

= 2.0 

for Z = 4.0 


r 

0 

0.039 

0.545 

0.094 

0.751 

0.158 

0.771 

1 

30 

.146 

.599 

.304 

.635 

.474 

.542 

60 

.294 

.713 

.601 

.762 

.946 

.629 


90 

.376 

.721 

.744 

.770 

1.179 

.654 


r 

0 

0.115 

0.663 

0.233 

0.820 

0.349^ 

0.788 

J 

30 

.176 

.694 

.371 

.792 

.583 

.695 


60 

.275 

.799 

.583 

.860 

.936 

.704 


190 

.316 

.828 

.661 

.869 

1.609 

.697 


r 0 

0.202 

0.750 

0.415 

0.922 

0.659 

0.851 

< 

30 

.236 

.730 

.479 

.866 

.764 

.771 

60 

.321 

.635 

.612 

.695 

.938 

.586 


[90 

.375 

.553 

.681 

.570 

1.008 

.470 
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TABLE I.- SCATTERED AND DISSIPATED POWER RATIOS - Continued 


(b) h = 4.0 


^0 

01, 

SR 

DR 

SR 

DR 

SR 

DR 

deg 

for 2 

= 1.0 

for 5 

5 = 2.0 

for Z 

= 4.0 



0 

0.050 

0.585 

0.122 

0.860 

0.218 

1.03 


) 

30 

.161 

.647 

.375 

.797 

.617 

.763 

0.2 


60 

.312 

.762 

.654 

.827 

1.069 

.697 



. 90 

.365 

.803 

.740 

.824 

1.154 

.653 



0 

0.138 

0.693 

0.298 

0.927 

0.483 

0.956 



30 

.202 

.703 

.437 

.872 

.742 

.828 

0.4 

< 

60 

.342 

.655 

.666 

.707 

1.038 

.606 



l90 

.414 

.597 

.765 

.604 

1.139 

.498 



r 0 

0.255 

0.639 

0.505 

0.803 

0.808 

0.801 

0.8 

} 

30 

.289 

.619 

.565 

.728 

.892 

.694 

< 

60 

.324 

.674 

.642 

.744 

1.005 

.652 



90 

.392 

.543 

.718 

.569 

1.069 

.499 


(c) h = 6.0 


1 


SR 

DR 

SR 

DR 

SR 

DR 

deg 

for Z 

= 1.0 

for Z 

^ = 2.0 

for Z 

= 4.0 



r ° 

0.065 

0.581 

0.161 

0.846 

0.288 

0.984 

0.2 


30 

.165 

.670 

.395 

.852 

.732 

.835 

\ 

60 

.320 

.758 

.666 

.815 

1.072 

.674 



L 90 

.390 

.747 

.760 

.773 

1.165 

,621 



0 

0.181 

0.590 

0.377 

0.728 

0.599 

0.697 

0.4 


30 

.236 

.628 

.496 

.754 

.828 

.717 

\ 

60 

.358 

.637 

.695 

.685 

1,076 

.592 



90 

.376 

.722 

.740 

.742 

1.137 

.596 



r 0 

0.272 

0.607 

0.540 

0.738 

0.865 

0.743 

0.8 


30 

.291 

.632 

.584 

.734 

.932 

.700 

j 

60 

.340 

.637 

.661 

.713 

1.031 

.649 



90 

V* 

.377 

.599 

.706 

.661 

1.074 

.605 
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TABLE I. - SCATTERED AND DISSIPATED POWER RATIOS - Concluded 


(d) h = 8.0 



Of, 

SR 

DR 

SR 

DR 

SR 

DR 

deg 

for Z 

= 1.0 

for 2 

= 2.0 

for 2 

; = 4.0 



0 

0.0692 

0.454 

0.171 

0.636 

0.308 

0.710 

0.2 


30 

.174 

.663 

.415 

.842 

.773 

.829 

\ 

60 

.359 

.667 

.708 

.716 

1.109 

.611 



90 

.440 

.609 

.814 

.596 

1.207 

.469 



r 0 

0.122 

0.509 

0.267 

0.704 

0.467 

0.773 

0.4 


30 

.261 

.596 

.531 

.713 

.869 

.704 

< 

60 

.343 

.689 

.682 

.748 

1.072 

.638 



90 

V 

.424 

.591 

.785 

.585 

1.165 

.466 



r ° 

0.188 

0.546 

0.411 

0.700 

0.726 

0.738 

0.8 

/ 

30 

.301 

.596 

.594 

.692 

.945 

.671 


60 

.337 

.637 

.660 

.713 

1.033 

.648 



L 90 

.340 

.658 

.668 

.724 

1.042 

.637 




(e) 

h = 10.0 






deg 

SR 

DR 

SR 

DR 

SR 

DR 


for Z 

= 1.0 

for 2 

= 2.0 

for Z = 4.0 



r " 

0.085 

0.554 

0.200 

0.807 

0.354 

0.978 

0.2 


30 

.186 

.631 

.432 

.795 

.790 

.795 

> 

60 

.384 

.595 

.737 

.615 

1.131 

.520 



L 90 

.429 

.648 

.800 

.664 

1.200 

.548 



r " 

0.193 

0.582 

0.414 

0.735 

0.688 

0.756 

0.4 

< 

30 

.256 

.587 

.533 

.694 

.885 

.656 

60 

.374 

.601 

.721 

.630 

1.106 

.530 



L 90 

.381 

.717 

.748 

.739 

1.150 

.596 



f 0 

0.289 

0.563 

0.578 

0.650 

0.920 

0.634 

0.8 


30 

.304 

.601 

.608 

.686 

.971 

.648 

s 

60 

.345 

.635 

.680 

.689 

1.064 

.608 



[90 

.338 

.716 

.687 

.766 

1.087 

.647 
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The ratio of scattered power to incident power (SR = varied from about 0.04 to 

\ ^ 1 / 

over 1.0. It was a minimum for grazing incidence (a = 0°) and a maximum for normal 
incidence {a = 90°) and increased monotonically with real impedance (resistance). As the 
resistance becomes infinite, SR must approach the value for the hard cylinder, A com- 
putation was made for the case of the hard cylinder and checked by using tabulated values 
of the radial functions (ref. 17). The numerical result for scattered power was found to 
be correct for h = 2,0 and a = 0°. 


The ratio of the dissipated power to the incident power ^DR = must have a max- 
imum value between Z = 0 and Z = «>. The value of DR depends on the product of 
the pressure drop and velocity at the cylinder surface. For zero resistance the pressure 
drop is zero, and for infinite resistance the velocity is zero whereas the pressure differ- 
ential must remain finite. It may be seen from the table that the optimum (maximum - 
dissipation) resistance is usually aroimd 2, except for grazir^ incidence on a thin = 0.2) 
cylinder where it is near 4. This large value of optimum resistance is a surprising result, 
since in one -dimensional problems, a resistive impedance of 1 gives the maximum absorp- 
tion, For thin cylinders DR was an increasing function of the angle of incidence for 
Z = 1.0 and a decreasing function of the angle for Z = 4.0. Another interesting result 
is that the maximum value for DR is just slightly larger than 1 for grazing incidence on 
a thin cylinder. Thus, it may be concluded that the maximum percentage of energy which 
can be removed from soimd propagating in a duct is equal to the percentage of duct area 
blocked by the splitter rings and spokes. 


In the case of normal incidence, the dissipation by an elliptic cylinder with eccen- 
tricity of nearly 1 should be roughly comparable to the dissipation by infinite parallel 
planes of porous material. If the length of the minor axis is held constant while the major 
axis is made much longer, then two sides of the cylinder will become nearly straight and 
parallel. This comparison is illustrated in figure 2. The data in this figure are for 

= 0.2, which corresponds to a cylinder whose major axis is roughly five times the length 
of the minor axis. The ratio of the dissipated power to the power in the normally incident 
plane wave is plotted against h, which may be interpreted as a frequency parameter, for 
the condition of normal Incidence. The curves in figure 2 give this ratio for parallel 
planes whose spacing is equal to the length of the minor axis of the cylinder. The maxi- 
mum dissipation and minimum dissipation for the parallel planes occur at very nearly the 
same frequencies as the maximum and minimum points for the cylinder, and the maxima 
have nearly the same values. 

Data for grazing incidence are also shown in figure 2. Only about one -fifth as much 
energy is dissipated from the grazing waves as from the normally incident wave. Dissi- 
pation ratios for the grazing incidence waves are shown in figure 3. The maximum and 


ji. 
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Z = 1.0 
Z = 2.0 
Z = 4.0 


Elliptic cylinder 


— Z = 1.0 

— Z = 2. 0 


Parallel planes 


= 0-2 


Normal 
(0-= 90°) 


-i 

2 sinh 4 

T 


Grazing 
(or= 0°) 


Figure 2.- Comparison of dissipation ty elliptic cylinder and parallel planes. 
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Figure 5*- Dissipation ratio for grazing Incidence. 
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minimum values of DR occur at frequencies which correspond to the frequencies at 
which the maxima and minima occur for the normal -incidence condition, a result which 
indicates that the minor -axis length is the critical dimension for grazing -wave energy 
dissipation. The maximum dissipation occurs when the minor axis is about one -fourth 
wavelength. 


CONCLUDING REMARKS 

The exact solution for scattering from a porous elliptic cylindrical shell with linear 
acoustic impedance has been found in terms of Mathieu fimction series. Coefficients of 
these series are solutions of four infinite matrix equations. Elements of the matrices 
depend on integrals of the periodic Mathieu functions and on the radial Mathieu fxmctions. 
Series expressions for the integrals have been fovmd. The Mathieu functions may be eval- 
uated in terms of known trigonometric and Bessel function series. By use of highly accu- 
rate modern computers, it has been possible to evaluate these functions for a range of 
parameters for which no tables are available. Thus, it has been possible to present results 
for a range of frequencies for which only formal solutions exist in previous literature. 

For a special case in which the acoustic impedance is a certain function of position 
on the cylinder, the infinite matrices reduce to diagonal forms, so that simple closed-form 
expressions in the series solutions are available. 

It is possible to obtain numerical results inside and outside of the cylinder in both 
the near field and the far field. Near afield solutions have been used to obtain the power 
dissipated in the porous shell. 

The periodic Mathieu functions used in this paper could easily be evaluated to eight 
significant figures. The radial Mathieu fimctions were more difficult to evaluate. The 
Bessel function series was superior for evaluating the radial Mathieu function of the first 
kind. This series had the property that higher order functions were more accurate for 
larger values of the parameter s (the square of the dimensionless focal length). Since 
s is proportional to frequency -squared, this property makes it possible to find solutions 
at the higher frequencies when the higher order functions are needed. 

The Bessel fimction product series was superior for computing the radial Mathieu 
function of the second kind. The accuracy of this series decreased with increasing s, or 
frequency. The complex radial Mathieu function was computed to six significant figures 
in the range 4 ^ s = 100. 

The integrals which give the coupling terms in the scattering equation were evaluated 
by double infinite series to eight significant figures. 

A study of scattering and dissipation of energy by the porous cylindrical shell has 
been made. Scattering is negligible for grazing incidence on the cylinder and is always 
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a maximum for normal incidence. The scattering increases monotonically with resistance 
of the porous material. 

The energy dissipated varies with frequency and resistance. The maximum dissi- 
pation usually occurs near a resistance of about 2; however, for grazing incidence the 
maximum dissipation occurs near a resistance of 4. 

The ratio of dissipated energy to incident energy reaches a maximum of about 1 for 
the grazing -incidence condition. This maximmn dissipation occurs when the minor axis 
of the cylinder is about one-fourth wavelength. This result suggests that the maximum 
percentage of energy which can be removed by placing a splitter in a duct is numerically 
equal to the percentage of duct area blocked by the splitter. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 21, 1971. 
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APPENDIX 


EVALUATION OF COUPLING COEFFICIENTS 

The integrals ^ and N^ from equations (23) to (26) are defined as 

p2TT Sei(v)Ser(v)dv , , 

r = ? r (38) 

’ ^ *^0 ^2(cosh 2 |q - cos 2r)) 

and 

.. _ 2 Soi(r})SOr(v)dv 

N? r ~ ^ \ - (39) 

^ ->0 |/2(cosh 2 ^q - cos 2rj) 

In order to evaluate the integrals, each term in the integrand wUl be expanded as a Fourier 
series in the variable t] as follows: 


00 


Se2r(??) = ^ 

(40) 

k=0 


00 

Se2r+l(^) = ^ cos(2k + l)v 

(41) 

k=0 


00 

®°2r+l("?) = Z °°2k+'l^^ ®^"(2k + 1 ) 7 ? 

(42) 

k=0 


CO 

So 2 r+ 2 (’ 7 ) = ^ °° 2 k+?^ + 2)77 

(43) 


k=0 


When these series are substituted into equations (38) and (39) for term-by-term integra- 
tion, it will be necessary to evaluate the integral 


I 


m 



cos iriT] drj 

|/2(cosh 2 |q - cos 277 ) 


(44) 
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APPENDIX - Continued 


The radical in the integrand of equation (44) may be represented as a series of 
Gegenbauer polynomials. Rainville (ref. 19) gives the generating function definition of 
these polynomials as 

00 

(l - 2xt +t2)"^= ^ C^(x)t’^ (45) 

n=0 


The denominator of the integrand in equation (44) may be written as 
|2(cosh - cos 2ri^ ^ = je^^°(l - 2 cos 2rje + e 


- 1/2 


Therefore 


(46) 


2(cosh 2 |q - cos 2 ? 7 jj = e ^ ^n^^ ^^^o 


n=0 


(47) 


The Gegenbauer polynomial in equation (47) may be given as (ref. 19, p. 283) 


n (I) (^] cos 2(n - 2k)Tj 


k=0 


(48) 


(l/'(l)(l) • • • (“-i)- 

Equations (47) and (48) combined give the desired expansion for the radical: 

cosh - cos 2^)] ^ = e-«° f f 

n=0 k=0 


[ 2 (, 


k;(n - k): 


(49) 


From equation (49) it may be seen that the integral in equation (44) is zero for all 
odd values of m; consequently, the only integrals which it is necessary to evaluate are 
of the form 


I 


2m “ 



cos 2m77 dri 
^2 (cosh 2 |q - cos 277 ) 


(50) 
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APPENDIX - Continued 


Substituting equation (49) into equation (50) gives 




^-2n5o y 2(n - 2k)?; cos 2mt) dtj 

n=m ktb - «•' 


( 51 ) 


That is, 


^ -(2m+l)4o Y (2)k(2)(m+2p-k) -4p?o o/ o o, ^ o ^ 

l2m = e TJT m TJ p - k lT ^ Jq 


(52) 


In equation (52) nonzero terms occur only when k = p or k = p+m; therefore. 


hm 


2ue 


/1\ /1\ 

(2m+l)|^ y \2)p\2)p^^ ^-4pgo 

L p:(p+m): 
p=0 


(53) 


Since the factorial terms in equation (53) form a decreasing sequence, the error 
due to truncating the series is 


E < 


pca 

Jn 


i\ /i\ 

p\2/p^m 


e ^^^Odp 


i\ /n 

2 /N\ 2 /N+m e' 


•4N^ 


pj(p + m): 


n;(n + m): Ho 


-4N^ 
e 


4^0 


(54) 


where N is the number of required terms. 

If an error less than e“^ is desired when 4^^ « 1, an estimate of the number of 
required terms is 


N > 


F - In Ho 

4^0 


(55) 


where F is the number of significant figures. With N determined by equation (55), 
equation (53) becomes 


i2m 


~ 27re 


N I-) 

■( 2 m+l)|o y U/pU/p+m ^-4pg( 

Zy p:(p+m): 
p=0 


(56) 
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APPENDIX — Continued 


Once the integrals 12m evaluated, they may be used to evaluate the integrals 
in equations (38) and (39). First, consider the product of the Mathieu functions in the 
integrand of equation (38). H Se^ and Se^. both have period tt, their product may be 
written as 


Se2^Se2r - ^ 
k=0 



cos 2k?7 cos 2jrj 


(57) 


Since cos 2kr7 cos 2j?7 = i cos 2(k - j)rj + ^ cos 2(k + j)p, the argument of each term in the 
series is an even multiple of p. Consequently, 


M, 


21, 2r h Z Z 
k=0 j=0 


^®2k 


j)g(2r ) 


^2|k-j I + ^2(k+j) 


(58) 


The product of Mathieu functions with different periods gives a series whose terms are 
all odd multiples of rj. It was noted previously that the integral of these terms is zero; 
therefore. 


'^21+1, 2r = 0 


(59) 


and 


'^21,2t+1 = 0 


(60) 


The product of even Mathieu functions with period 2tt is 


®®2Z+l®®2r+l - 2 
k=0 


oo 




De 


(2^+l) 

2k+l 


De 


(2r+l) 

2j+l 


cos(2k + 1)77 cos(2j + 1)77 


(61) 


Therefore 


M. 


2Z+l,2r+l “ h 

k=0 j=0 


De 


(2Z+1) 

2k+l 


De 


(2r+l) 

2j+l 


^2lk-ij 


+ ^2(k+j+l) 
2 


(62) 
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APPENDIX — Concluded 


A similar evaluation for integrals involving odd Mathieu fimctions gives 

2 V V r^ (2^+1) (2r+l) 

N2Z+l,2r+l = h 2/ Z ^°2k+l ^°2j+l 

k=0 j=0 

OO 

TM 2 V V n^(2Z+2) ^ (2r+2) 

N2Z+2,2r+2 = h 2/ Z ^ 2k+2 °°2j+2 

k=0 j=0 


^2|k-j| - ^2(k+j+l) 
2 

^2|k-j| - ^2(k+j+2) 
2 
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SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. S0546 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs^ 
Technology Utilization Reports and 
Technology Surveys. 


